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The measured ft-vahies for superallowed 0^ — > 0^ nuclear /3-decay can be used to obtain the 
value of the vector coupling constant and thus to test the unitarity of the Cabibbo-Kobayashi- 
Maskawa matrix. An essential requirement for this test is accurate calculations for the radiative 
and isospin symmetry-breaking corrections that must be applied to the experimental data. We 
present a new and consistent set of calculations for the nuclear-structure-dependent components of 
these corrections. These new results do not alter the current status of the unitarity test - it still 
fails by more than two standard deviations - but they provide calculated corrections for eleven new 
superallowed transitions that are likely to become accessible to precise measurements in the future. 
The reliability of all calculated corrections is explored and an experimental method indicated by 
which the structure-dependent corrections can be tested and, if necessary, improved. 

PACS numbers; 23.40.Bw, 23.40.Hc 



I. INTRODUCTION 

Superallowed 0+ 0+ nuclear /3-decay depends 
uniquely on the vector part of the weak interaction. 
When it occurs between T — 1 analog states, a precise 
measurement of the transition /i-value can be used to 
determine Gv, the vector coupling constant. This result, 
in turn, yields Vud, the up-down element of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. At this time, it is 
the key ingredient in one of the most exacting tests avail- 
able of the unitarity of the CKM matrix, a fundamental 
pillar of the minimal Standard Model. 

Currently, there is a substantial body of precise ft- 
values determined for such transitions and the exper- 
imental results are robust, most input data having 
been obtained from several independent and consistent 
measurements ^. In all, /i- values have been de- 
termined for nine 0+ 0+ transitions to a preci- 
sion of ~ 0.1% or better. The decay parents - ^°C, 
"O,26™Al,34cl,38™K,42Sc,46v,50Mn and ^^Go - span a 
wide range of nuclear masses; nevertheless, as anticipated 
by the Conserved Vector Current hypothesis, CVC, all 
nine yield consistent values for Gy, from which a value 
of 



Vud = 0.9740 ± 0.0005 



(1) 



is derived. The unitarity test of the CKM matrix, made 
possible by this precise value of V^d, fails by more than 
two standard deviations [Q: viz. 

yld + yls + = 0.9968 ± o.oow. (2) 

In obtaining this result, we have used the Particle Data 
Group's B recommended values for the much smaller 



'Present address: Department of Physics, Queen's University, 
Kingston, Ontario K7L 3N6, Canada 



matrix elements, Vus and Vuh- Although this deviation 
from unitarity is not completely definitive statistically, 
it is also supported by recent, less precise results from 
neutron decay If the precision of this test can be 
improved and it continues to indicate non-unitarity, then 
the consequences for the Standard Model would be far- 
reaching. 

The potential impact of definitive non-unitarity has led 
to considerable recent activity, both experimental and 
theoretical, in the study of superallowed 0+ 0+ tran- 
sitions, with special attention being focussed on the small 
correction terms that must be applied to the experimen- 
tal /i- values in order to extract Gv. Specifically, Gv is 
obtained from each /t-value via the relationship ^ 



K 



2G2(1 + Av)^ 



(3) 



with 



= (8120.271 ± 0.012) X 10"^"GeV"*s, (4) 

where / is the statistical rate function, t is the partial 
half-life for the transition, 5c is the isospin-symmetry- 
breaking correction, 5r is the transition-dependent part 
of the radiative correction and is the transition- 
independent part. Here we have also defined J-t as the 
"corrected" /i-value. 

It is now convenient to separate the radiative correc- 
tion into two terms: 



5r=^ 5'r, + Sns 



(5) 



where the first term, (5^, is a function of the electron's en- 
ergy and the charge of the daughter nucleus, Z; it there- 
fore depends on the particular nuclear decay, but is in- 
dependent of nuclear structu re. T he second term, Sns, is 
discussed more fully in Sec. [IB but its evaluation, like 



2 



6c, depends on the details of nuclear structure. To em- 
phasize the different sensitivities of the correction terms 
we re- write the expression for Tt as 



Tt = ftil + d'ji)il + 6Ns-Sc) (6) 



where the first correction in brackets is independent 
of nuclear structure, while the second incorporates the 
structure-dependent terms. The term S'j^ has been calcu- 
lated from standard QED, and is currently evaluated to 
order Za^ and estimated in or dcr Z^a^ [|, |; its values 
are around 1.4% and can be considered very reliable. The 
structure-dependent terms, Sns and 6c, have also been 
calculated in the past but at various times over three 
decades and with a variety of different nuclear models. 
Their uncertainties are larger. This paper specifically 
addresses these correction terms with a view to reducing 
their uncertainties. 

Though depending on the nuclear shell-model, calcu- 
lations for (5jvs and Sc have been carefully linked to 
other related observables such as the neutron and proton 
binding energies, the b- and c-coefficients in the Isobaric 
Multiplet Mass Equation (IMME), and the non-analog 
0+ 0+ transition rates (see, for example, ^). 
Given this linking to observables and the more general 
success of the shell model in this mass region, calculations 
of Sns and Sc should also be rather reliable. Neverthe- 
less, conservative uncertainties have been applied - they 
are of order 0.1% {i.e. ^ 10% of their own value) - and 
these become major contributors to the overall uncer- 
tainty on the unitarity test. To illustrate: the uncertainty 
obtained for Vud in Eq.(|l])is ±0.0005; the contributions 
to this uncertainty are 0.0001 from experiment, 0.0001 
from S'l^, 0.0003 from Sc - Sns, and 0.0004 from A^. 
If the unitarity test is to be sharpened, then the most 
pressing objective must be to reduce the uncertainties on 
and {Sc — Sns)- The latter is clearly the most im- 
portant area where nuclear physics can play an critical 
role. There is considerable activity, both experimental 
and theoretical, now underway in probing these nuclear- 
structure-dependent corrections with a view to reducing 
the uncertainty that they introduce into the unitarity 
test. 

Since the goal of experiments will generally be to test 
and constrain the calculated structure-dependent correc- 
tions, an important first step is to have a set of consistent 
calculations that apply both to the nine well-known tran- 
sitions already used for the unitarity test and to possi- 
ble new cases yet to be studied. In what follows, we 
present new calculations of Sc and Sns, in which con- 
sistent model spaces and approximations have been used 
for both correction terms and for a large repertoire of su- 
perallowed transitions, new and old. These will provide a 
consistent standard for future experimental comparison. 



II. THEORETICAL CORRECTIONS TO 
SUPERALLOWED DECAYS 

As described in the introduction, there are four the- 
oretical correction terms involved in extracting Vud 
from experimental /i-values: the radiative corrections 
that are independent of nuclear structure ((5^ and 
A^), the nuclear-structure-dependent radiative correc- 
tion {Sns) and the isospin-symmetry-breaking correction 
{Sc). Though we will later present new calculations of 
the last two, in this section we present an overview of all 
four terms. This overview is placed in the context of a 
unitarity test that has failed by more than two standard 
deviations. In particular, we assess whether the failure to 
meet unitarity can be removed by plausible adjustments 
in these calculated corrections. What changes would it 
take to restore unitarity? For example, S'j^, would have 
to be shifted downwards by 0.3% {i.e. as much as one- 
quarter of its current value) for all nine currently well- 
measured nuclear transitions; or {Sc — Sns) would have 
to be shifted upwards by 0.3% (over one-half their value), 
for all nine cases; or some combination of the two. We 
will argue that such shifts are very improbable. 

A. Radiative corrections independent of nuclear 
structure 

The radiative correction comprises a transition- 
dependent term, Sr, and a transition-independent term, 
A^. The transition-dependent term is further divided 
into S'ji, which does not depend on nuclear structure, and 
Sns, which is structure dependent. We consider first the 
structure-independent terms, which are written: 

S'r = ^[g{E^) + S, + Ss], (7) 
zvr 
a 

Al = — [41n(TOz/mp) + ln{mp/nij,) + 2CBorn] H , 

ZTT 

where the ellipses represent further small terms of order 
0.1%. In these equations, E^^ is the maximum electron 
energy in beta decay, ruz the Z-boson mass, nip the pro- 
ton mass, ttia the oi-meson mass, and S2 and ^3 the order- 
Za^ and -Z'^c? contributions respectively. The function 
g{Ee, Ern), which depends on the electron energy, was 
first defined by Sirlin pO| as part of the order-a univer- 
sal photonic contribution arising from the weak vector 
current; it is here averaged over the electron spectrum 
to give 'g{Em). Finally, the term Cboiii comes from the 
order-a axial- vector photonic contributions. 

Calculated values for all three components of (5^ are 
given in Table ^ There have been two independent cal- 
culations 1^, ^ of both S2 and Sz; they are completely 
consistent with one another if proper account is taken 
of finite-size effects in the nuclear charge distribution. 
The values listed in Table p|are our recalculations |^ us- 
ing the formulas of Sirlin B but incorporating a Fermi 
charge-density distribution for the nucleus. Note that we 
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TABLE L Calculated nucleus-dependent radiative correction, 
Sji, in percent units, and the component contributions as iden- 
tified in Eq. (|). 











Or 


T, = -1: 












1.468 


0.182 


0.005 


1.65(1) 


14 Q 


1.286 


0.227 


0.008 


1.52(1) 


^**Ne 


1.204 


0.268 


0.013 


1.48(1) 


99 n r 


1.121 


0.305 


0.018 


1.44(2) 


^''Si 


1.055 


0.338 


0.024 


1.42(2) 


30 g 


1.005 


0.363 


0.030 


1.40(3) 


''^Ar 


0.963 


0.392 


0.037 


1.39(4) 


''^^Ca 


0.928 


0.417 


0.044 


1.39(4) 


d.9 m. 

42 rp^ 


0.906 


0.449 


0.053 


1.41(5) 


= 0: 










26m ^1 


1.110 


0.325 


0.021 


1.46(2) 


34 (^j 


1.002 


0.388 


0.034 


1.42(3) 


38m 


0.964 


0.413 


0.041 


1.42(4) 


*2Sc 


0.939 


0.448 


0.049 


1.44(5) 


46y 


0.903 


0.468 


0.057 


1.43(6) 


5»Mn 


0.873 


0.494 


0.065 


1.43(7) 


54 Co 


0.843 


0.507 


0.073 


1.42(7) 


^^Ga 


0.805 


0.567 


0.091 


1.46(9) 


'^''As 


0.791 


0.589 


0.100 


1.48(10) 


^«Br 


0.777 


0.609 


0.110 


1.50(11) 


^*Rb 


0.763 


0.627 


0.120 


1.51(12) 



have followed Sirlin in assigning an uncertainty equal to 
{a/2Tr)d3 as an estimate of the error made in stopping 
the calculation at that order. 

To assess the changes in S'j^ that would be required in 
order to restore unitarity, it is helpful to rewrite Eq. (Q) 
in terms of the typical values taken by its components: 
viz. 

(5^- 1.00 + 0.40 + 0.05%, (8) 

If the failure to obtain unitarity in the CKM matrix with 
Viid from nuclear beta decay is due to the value of this 
term alone, then S'j^ must be reduced to 1.1%. This is not 
likely. The leading term, 1.00%, involves standard QED 
and is well verified. The order-Za^ term, 0.40%, while 
less secure has been calculated twice ^ independently, 
with results in accord. 

Taking a similar approach for the nucleus-independent 
radiative correction, we write 

2.12 - 0.03 + 0.20 + 0.1% ~ 2.4%, (9) 

of which the first term, the leading logarithm, is unam- 
biguous. Again, to achieve unitarity of the CKM matrix, 
would have to be reduced to 2.1%: i.e. all terms 
other than the leading logarithm must sum to zero. This 
also seems unlikely. The adopted value of the nucleus- 
independent radiative correction has been set at |l2] 



A;^ = (2.40 + 0.08)%. (10) 

Note this value differs slightly (but within errors) from 
an earlier value ||ll[ because of the decision by Sirlin |lj] 
to centre the cut-off parameter rrij^, where {niai/'i.) < 
TOa £ 2mai J exactly at the ai-meson mass when eval- 
uating the axial contribution to the radiative-correction 
loop graph. This range of possible values for ttia is the 
dominant contributor to the error in Eq. ([l0|). 

B. The 5ns correction 

The nuclear-structure-dependent part of the radiative 
correction is denoted 5ns- Although for the superallowed 
transition we are discussing a purely vector interaction 
between spin 0+ states, the axial-vector interaction does 
play a role in the radiative corrections. An axial-vector 
interaction may flip a nucleon spin and then be followed 
by an electromagnetic interaction that flips it back again. 
This axial contribution, denoted C, can be further di- 
vided into two terms depending whether the weak and 
electromagnetic interactions occur on the same nucleon 
or on two separate nucleons: 

C = Ceorn + Cns, 
5ns = —Cns- (11) 

TT 

Here Cboih refers to the Born graph in which the axial- 
vector and electromagnetic interactions occur on the 
same nucleon. This term is universal - i. e. the same in all 
nuclei - so it is not included in Sns but is placed in the 
nucleus-independent radiative correction A^, see Eq. (|^). 
The term, Cns 7 refers to the case when the axial- vector 
and electromagnetic interactions occur on different nucle- 
ons. The calculation of this term depends on the details 
of nuclear structure. 

In the earliest calculations of 5ns |l^ , the axial- 
vector and electromagnetic vertices were evaluated with 
free- nucleon coupling constants. Yet there is ample evi- 
dence in nuclear physics that coupling constants for spin- 
flip processes are quenched in the nuclear medium. Sub- 
sequently, Towner [|l6j revised his earlier results |jl^ us- 
ing quenching factors that had been obtained previously 
p7| , p^ , from studies of weak and electromagnetic 
transitions in nuclei throughout the region 10 < ^ < 54. 
These quenching factors depend weakly on both mass 
and shell-model orbital. 

There is a further consideration. The presence of 
quenching also breaks the universality of the Born term, 
C'sorn- Writing the evaluation of CBom with free- nucleon 
coupling constants as CBom (free) , then CBom (quenched) 
can be written: 

Cboiii (quenched) = 9CBom(free) (12) 
= CBom(free) + (q - l)CBom(free), 
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TABLE 11: Shell-model calculations of the nuclear-structure dependent component of the radiative correction, Sns- The four 
components that are summed to give Cns characterize the four electromagnetic couplings: os = orbital isoscalar, ss — spin 
isoscalar, ov — orbital isovector, and sv — spin isovector. 



Parent 


Unquenched 




Quenched Cns 


i1- 


l)x 


5ns {%) 


nucleus 


Cns 


OS 


ss ov sv 


total Csorn 


(free) 


Quenched Adopted 



= -1: 



10^ 


-1.669 


0.002 


-0.283 


-0.002 


-1.065 


-1.348 


-0.188 


-0.357 


-0.360(35) 


14 Q 


-1.360 


-0.008 


-0.341 


0.082 


-0.782 


-1.049 


-0.221 


-0.295 


-0.250(50) 


i»Ne 


-1.531 


-0.011 


-0.249 


-0.119 


-0.812 


-1.191 


-0.210 


-0.325 


-0.290(35) 


22 Mg 


-1.046 


-0.009 


-0.222 


-0.067 


-0.497 


-0.796 


-0.226 


-0.237 


-0.240(20) 


26Si 


-0.986 


-0.007 


-0.224 


-0.086 


-0.424 


-0.741 


-0.242 


-0.228 


-0.230(20) 


30 g 


-0.800 


0.002 


-0.287 


0.020 


-0.300 


-0.566 


-0.257 


-0.191 


-0.190(15) 


3*Ar 


-0.770 


0.014 


-0.322 


0.061 


-0.272 


-0.519 


-0.273 


-0.184 


-0.185(15) 


3«Ca 


-0.693 


0.041 


-0.358 


0.091 


-0.214 


-0.440 


-0.288 


-0.169 


-0.180(15) 


42^i 


-1.011 


-0.016 


-0.181 


-0.225 


-0.354 


-0.776 


-0.256 


-0.240 


-0.240(20) 



r, =0: 



26m ^1 


0.352 


-0.007 


-0.224 


0.086 


0.424 


0.279 


-0.242 


0.009 


0.009(20) 


34C1 


-0.135 


0.015 


-0.333 


-0.064 


0.280 


-0.101 


-0.273 


-0.087 


-0.085(15) 


38m 


-0.276 


0.042 


-0.363 


-0.093 


0.216 


-0.198 


-0.288 


-0.113 


-0.100(15) 


*2Sc 


0.472 


-0.016 


-0.182 


0.228 


0.358 


0.389 


-0.256 


0.031 


0.030(20) 


46y 


0.101 


-0.004 


-0.197 


0.099 


0.198 


0.096 


-0.263 


-0.039 


-0.040(7) 


^''Mn 


0.054 


-0.009 


-0.184 


0.104 


0.152 


0.063 


-0.270 


-0.048 


-0.042(7) 


54 Co 


0.161 


-0.013 


-0.180 


0.133 


0.203 


0.144 


-0.277 


-0.031 


-0.029(7) 


«2Ga 


0.172 


0.005 


-0.289 


-0.058 


0.445 


0.103 


-0.289 


-0.043 


-0.040(20) 


'^''As 


0.124 


0.006 


-0.291 


-0.070 


0.421 


0.066 


-0.295 


-0.053 


-0.050(20) 


™Br 


0.077 


0.009 


-0.295 


-0.083 


0.401 


0.032 


-0.301 


-0.063 


-0.060(20) 


^^Rb 


0.155 


0.009 


-0.261 


0.006 


0.353 


0.106 


-0.306 


-0.046 


-0.065(20) 



where q is the factor by which the product of the weak 
and electromagnetic couphng constants is reduced in the 
medium relative to its free-nucleon value. The first term 
in Eq. ( p^ remains universal, while the second term is 
now part of the nuclear-structure dependence of the ra- 
diative correction. Thus Sns is written 



Sns = - {Cats (quenched) + {q- l)CBorn(free)} . (13) 

TT 

We have calculated the S ns correction for a wide range 
of nuclei with 0^ (T = 1) ground or isomeric states 
that decay by superallowed /3-emission; we used the shell 
model with effective interactions as described in Ap- 
pendix A. Results for both quenched and unquenched 
coupling constants are given in Table All but the last 
column in that table give the results from one particular 
calculation for each parent nuclide. (In most cases, two 
or three independent calculations were performed for a 
single parent, each with a different shell-model Hamilto- 
nian.) The last column lists the values we adopt for Sns- 
these values result from our assessment of the quenched 
results from all calculations made for each decay ~ not 
just the ones shown in the previous columns - with un- 
certainties chosen to encompass the spread in the results 
from those calculations. 



Extra details are also given in columns 3-6 of the ta- 
ble for the quenched calculation since this is the version 
that we ultimately use in evaluating Vud- With two-body 
operators there are two types of contributions: those in 
which both interacting nucleons are in the valence model 
space, and those in which one nucleon is in the valence 
space and one is in the closed-shells core. In the lat- 
ter case a sum is required over all the core nucleons. 
The isospin structure of the operator is interesting to 
note: the weak interaction contribution is isovector, while 
the electromagnetic contribution is isoscalar or isovector. 
The combined operator therefore is either isovector or 
isotensor. (An isoscalar combination is just proportional 
to the unit operator in isospin space and does not induce 
a Fermi transition.) Both the valence nucleons and those 
in the core contribute to the result for isovector opera- 
tors, only the valence nucleons contribute to the isotensor 
operators. 

In Table ^ we show contributions to Cns from the 
various components of the electromagnetic interaction: 
orbital isoscalar (os), spin isoscalar (ss), orbital isovector 
(ov) and spin isovector (sv). Note that the spin contri- 
butions are larger than the orbital contributions. Fur- 
ther, and more interesting, the isoscalar and isovector 
contributions are in phase when the decaying nucleus has 
Tz = —1 and out of phase when the decaying nucleus has 
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Tz = 0. This indicates that much larger corrections are 
obtained in the = — 1 series than in the = se- 
ries. If one looks at mirror transitions, this effect alone 
contributes between 0.1% to 0.3% to a mirror asymme- 
try in the /t-values. Since current experiments aim at 
0.1% accuracy, this effect might just be at the edge of 
detectability. 



Also, with charge-dependent terms in the effective inter- 
action, the Fermi matrix elements to other non-analogue 
0+ states in the daughter are no longer exactly zero. For 
example, there could be small (usually less than 0.1%) 
branches to those excited 0+ states that are energeti- 
cally accessible to beta-decay. For the first excited (non- 
analog) 0+ state, we can write 



C. Isospin symmetry- breaking corrections 

Turning, next, to the isospin-symmetry breaking cor- 
rection, 6c, it too can be separated into two components: 



(16) 



<5c = <5ci + Sc'2- 



(14) 



The first term, Sci, arises from Coulomb and charge- 
dependent nuclear interactions that induce configuration 
mixing among the 0"*" state wave functions in both the 
parent and daughter nuclei. Being charge dependent, this 
mixing serves to break isospin symmetry between the 
analog parent and daughter states of the superallowed 
transition. The second term, dc2, is due to small dif- 
ferences in the single-particle neutron and proton radial 
wave functions, which cause the radial overlap integral 
of the parent and daughter nucleus to be less than unity. 
Strictly speaking, these two aspects of the calculation of 
(5(7 cannot be separated, but in all but one calculation 
to date (including those reported here) this division has 
been made. (The exception is the (0-1-2-1- A)hLj large- 
basis shell-model calculation of Navratil et al. [2^] for the 
lightest superallowed emitter, ^°C.) This division is akin 
to the division made in setting up a shell-model calcu- 
lation, where the configuration space is divided into a 
small, tractible valence space and a remaining excluded 
space. Then Sci arises from the charge-dependent mixing 
within the valence space, while 6c2 represents the conse- 
quence of mixing between configurations in the valence 
space with those in the excluded space; this consequence 
being manifested by a change in the single-particle radial 
wave function of the valence nucleons. 



1. The 3ci correction 

If, in a shell-model calculation, the effective interac- 
tion is isospin invariant, then the wave functions for the 
parent and daughter analogue states are identical, and 
the square of the Fermi matrix element between them 
(for isospin T = 1 states) is exactly — 2. In 

addition, beta transitions to all other 0+ states in the 
daughter are strictly forbidden. However, the addition 
of charge-dependent terms to the effective interaction 
causes the breaking of analogue symmetry. Under these 
conditions, the Fermi matrix element departs slightly 
from its isospin-invariant value. We write 



In a model calculation in which there are only two basis 
states, the depletion of Fermi strength in the ground- 
state transition is entirely picked up by the transition to 
the excited non-analogue 0+ state. Thus, 



(17) 



Further, if only two-state mixing is considered, the mag- 
nitude of Sci is inversely proportional to the square of 
the excitation energy of the excited 0+ state: i.e. 



Sci oc 



(18) 



2(1- Sci). 



(15) 



For our calculations, in which a large number of basis 
states play a role, Eqs. ( p7| ) and (|8|) are no longer exact. 
Even so, they remain approximately true and continue 
to be a useful guide. 

Calculations of Sci turn out to be very sensitive to 
the details of the model calculation. This would be a 
very unfortunate property if we were not able to adopt 
certain strategies that act to reduce the model depen- 
dence considerably. Because of the variation of Sci with 
(Ai5)^ (see Eq. (|8|)), it is important that the isospin- 
independent Hamiltonian produce a good quality spec- 
trum of 0"*" states. Since this is not always possible 
to achieve in the shell model, especially for nuclei near 
to closed shells, our first strategy is to compensate for 
this by scaling the calculated Sci values by a factor 
(Ai?)jjjj,Q/(Ai?)^xp^, the ratio of the square of the ex- 
citation energy of the first excited 0+ state in the model 
calculation to that known experimentally. The second 
strategy we adopt to reduce the model dependence was 
first used by Ormand and Brown 0, |2^. We constrain 
the charge-dependent part of the effective interaction to 
reproduce other charge-dependent properties of the 0"*" 
states, namely the coefficients of the isobaric multiplet 
mass equation (IMME)Q. 

There are three ways in which charge dependence en- 
ters our shell- model calculation. First, the single-particle 
energies of the proton orbits are shifted relative to those 
of the neutrons. The amount of shift is determined from 
the spectrum of single-particle states in the closed-shell- 
plus-proton versus the closed-shell-plus-neutron nucleus, 
where the closed shell is taken to be the nucleus used as 
a closed-shell core in that particular shell-model calcula- 
tion. These single-particle shifts are taken from experi- 
ment and are not adjusted. Second, a two-body Coulomb 
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TABLE III: Shell-model calculations of the isospin symmetry-breaking correction, Sci ■ 



Parent 
nucleus 



Measured IMME coefficients |21| 



5cii%) 



expt 


SM 


1 1 n "Spa 1 pH 


ciPfi 1 pH 


A H n"n1"prl 


6.18 


11.05 


0.002 


0.007 


0.010(10) 


6.59 


6.64 


0.049 


0.050 


0.050(20) 


3.63 


3.80 


0.212 


0.232 


0.230(30) 


6.24 


6.34 


0.010 


0.010 


0.010(10) 


3.59 


4.96 


0.022 


0.042 


0.040(10) 


3.79 


3.86 


0.186 


0.193 


0.195(30) 


3.92 


3.91 


0.031 


0.030 


0.030(10) 


3.38 


3.21 


0.026 


0.023 


0.020(10) 


1.84 


3.60 


0.065 


0.249 


0.220(100) 



T, = -1: 

*o 

^Ne 
^Si 
*Ar 



10( 
14( 
18j 
22j 
26 c 
30 c 
34^ 

3«Ca 

42^. 



-1.546 

-2.493 

-3.045(1) 

-3.814(1) 

-4.535(2) 

-5.185(2) 

-5.777(2) 

-6.328(3) 

-6.712(3) 



0.362 

0.337 

0.347(1) 

0.315(1) 

0.302(2) 

0.275(2) 

0.286(2) 

0.284(3) 

0.287(3) 



= 0: 



26m ^1 


-4.535(2) 


0.302(2) 


3.59 


4.96 


0.022 


0.041 


0.040(10) 


34(^1 


-5.777(2) 


0.286(2) 


3.92 


3.91 


0.103 


0.103 


0.105(20) 


38m 


-6.328(3) 


0.284(3) 


3.38 


3.21 


0.099 


0.089 


0.100(20) 




-6.712(3) 


0.287(3) 


1.84 


3.60 


0.019 


0.072 


0.060(30) 


46y 


-7.327(10) 


0.276(11) 


2.61 


3.92 


0.043 


0.097 


0.095(20) 


5«Mn 


-7.892(30) 


0.259(30) 


3.69 


4.23 


0.048 


0.063 


0.055(20) 


s^Co 


-8.519(25) 


0.276(25) 


2.56 


2.26 


0.058 


0.045 


0.040(15) 


^^Ga 


-9.463(70) 


0.265(25)" 


2.33 


2.26 


0.350 


0.330 


0.330(40) 


««As 


-9.95(15) 


0.262(25)" 


2.17' 


1.81 


0.356 


0.247 


0.250(40) 


'°Br 


-10.48(23) 


0.260(25)" 


2.01 


1.72 


0.479 


0.352 


0.350(40) 




-10.82(25) 


0.258(25)" 


0.508 


0.523 


0.122 


0.129 


0.130(60) 



"Estimated values extrapolated from a fit to c coefficients in 0+ 
states in A = 4n + 2 nuclei, 10 < A < 58. 

''Estimated value taken to be an average of the excitation energies 
of 0+ states in ^^Zn and ''"Se. 



interaction is added among the valence protons. The 
strength of this interaction is adjusted so that the b- 
coefficient of the IMME is exactly reproduced. Third, we 
add a charge-dependent nuclear interaction by increasing 
all the T = 1 proton-neutron matrix elements by about 
2% relative to the neutron-neutron matrix elements. The 
precise amount of this increment was determined by re- 
quiring that the c-coefficient of the IMME be exactly 
reproduced. 

For each of the nuclei appearing in the previous ta- 
bles, we list in Table |l| the values of the corresponding 
measured IMME coefficients, b and c, together with the 
known excitation energy, £^^(0+), of the lowest excited 
0+ state in their daughters. As already explained, all 
our shell-model calculations were adjusted to reproduce 
exactly the values of b and c, and any discrepancy be- 
tween the calculated and experimental values of i?2;(0+) 
was compensated for by scaling the calculated results for 
Sci- As we did in Table ||, Columns 5-7 of in this table 
give the results from one particular calculation for each 
parent nucleus. These columns list the calculated 0^ ex- 
citation energy and Sci values, both unsealed and scaled 
for any E^lO'^) discrepancy. Finally, the eighth column 
gives the 5ci values we adopt. These values result from 



our assessment of the results of all calculations made for 
each decay - not just the ones shown in columns 5-7 - 
with uncertainties chosen to encompass the spread in the 
results from those calculations and to include the uncer- 
tainty in the IMME b- and c-coefhcients. 



For the nuclei with A > 38 there are excited (non- 
analogue) O"*" states in the daughter nuclei that are ac- 
cessible to beta decay. Some of the Fermi transitions to 
these states have also been measured [0, . In Table fv 



we list one set of calculated S^^ values, both unsealed 
and scaled, along with the value of Sq^ we adopt based 



on the same assessment as that described for Table III 



As before, the assigned errors reflect both the spread 
among the different calculations and the uncertainties in 
the IMME coefficients. The measured branching ratios 
were then converted to S^^ values (see Eq. (^6|)), which 
appear in the last column of the table. With the possible 
exception of the results for ^°Mn, the agreement between 
theory and experiment is entirely satisfactory. 
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TABLE IV: Shell-model calculations for the square of the 
Fermi matrix clement to the first excited 0"*^ state, Sqi. 



Vc(r) = Ze'^/r, 



Parent 
nucleus 






5hi{%) 




unsealed 


scaled 


adopted 


expt 












38m 


0.068 


0.062 


0.090 (30) 


< 0.28" 


^^Sc 


0.007 


0.029 


0.020 (20) 


0.040(9)' 


46y 


0.020 


0.046 


0.035 (15) 


0.053(5)" 


5«Mn 


0.038 


0.049 


0.045 (20) 


< 0.016" 


54 Co 


0.049 


0.038 


0.040 (20) 


0.035(5)" 


e^Ga 


0.089 


0.084 


0.085 (20) 




««As 


0.027 


0.019 


0.020 (20) 




^«Br 


0.095 


0.070 


0.070 (20) 






0.045 


0.047 


0.050 (30) 





"From Hagberg et al. [g| 
''From Daehnick and Rosa 



averaged with earlier results. 



2. The 5c2 correction 

The second isospin symmetry-breaking correction, Sc2 , 
accounts for the difference in radial forms between the 
proton in the parent /3-decaying nucleus and the neutron 
in the daughter nucleus. These radial forms are inte- 
grated together and, if there were no difference between 
them, the integral would just be the normalization inte- 
gral of value one. The departure of the square of this 
overlap integral from unity corresponds to 6c2- There is 
a strong constraint on any calculation of Sc2' the asymp- 
totic forms of the radial functions must be matched to 
the separation energies, Sp and 5„, where Sp is the pro- 
ton separation energy in the parent nucleus and Sji is 
the neutron separation energy in the daughter nucleus. 
These separation energies are well known and and may 
be found in any atomic mass table. It is the size of the 
difference between Sp and Sn and whether or not the ra- 
dial wave functions have nodes that principally determine 
the magnitude of 6c2- 

Our calculations of this correction follow closely the 
methods described in our earlier work We use a 
Saxon- Woods potential defined for a nucleus of mass A 
and charge Z + 1 as: 



V{r) = ~Vof{r)-V,gir)l.(T + Vc{r) 
-Vgg{r)-Vnh{r), 



(19) 



for r > Rc 

2 \ 



for r < Rr 



(20) 



with R = ro{A - l)!/^ and R^ = r,{A - ly/^. Note 
that g(r) is rendered dimensionless through the use of 
the pion Compton wavelength, {h/raT^c) = 2 fm^. The 
first three terms in Eq. (|l9) are the central, spin-orbit 
and Coulomb terms respectively. The fourth and fifth 
terms are additional surface terms whose role we discuss 
shortly. The parameters of the spin-orbit force were fixed 
at standard values, Vs — 7 MeV, Vg — 1.1 fm and = 
0.65 fm, leaving four parameters to be determined: i?c, 
the radius of the Coulomb potential, and Vb, ro and a 
characterizing the strength, range and diffuseness of the 
Saxon- Woods potential. 

To determine the radius of the Coulomb potential, Rc, 

we first obtained the charge mean-square radius, {'r'^)\(^ , 
of the decaying nucleus. We used results from electron 

which actually provide the 



24 



scattering experiments 
charge radius of a stable isotope of each element rather 
than the beta-decaying isotopes of interest here. How- 
ever, by examining the data on isotope shifts of charge 
radii we could make corrections for this effect to arrive 
at radius values applicable to the decaying nuclides; we 
enlarged the assigned error accordingly. Our selected 
values of {r'^)]l^ and their assigned errors are listed in 
Table To obtain an appropriate value for i?c, two fur- 
ther adjustments are required to the experimental values 
of {r'^)ll^- first, the finite size of the proton must be 
incorporated and second, because the shell model uses 
A single-particle coordinates rather than (^4 — 1) relative 
coordinates, a centre-of-mass correction must be applied. 
With a Gaussian form for the proton single-particle den- 
sity and harmonic oscillator wave functions for the shell 
model, the shell-model rms radius, {r'^)lm, relates to the 
experimentally measured rms radius via 



{r'U = (-^>.n. + Ual- b^A) , 



(21) 



where 



where = 0.694 fm is the length parameter in the pro- 
ton density and b is the length parameter of the har- 
monic oscillator, approximately b^ ~ A^l'^ fm-^. The 
Coulomb potential in Eq. (|2^) is that of a uniformly 
charged sphere. We match the charge radius of this dis- 
tribution with (r^)sm^ to determine the radius, i?c. 



/(r) = {l-|-exp((r-i?)/a)}-\ 

M (r^R, 
exp 



.9(0 



X ^ 1 + exp 



r-R. 



R^ = -(r^) 



(22) 



Finally, it remains to determine the parameters of the 
central potential, Vb/(r). The diffuseness is fixed at the 
same value as that of the spin-orbit potential, a = 0.65 
fm, for all A values except the lightest, A — W> and 14, 
for which we used a = 0.55 fm. The well depth, Vqj 
was adjusted case-by-case so that the asymptotic form 
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TABLE V: Calculations of 5c2 with Saxon- Woods radial functions, without parentage expansions, 5^2^ ^^nd with parentage 
expansions, 5q2i ^ci and Sc^. 



Parent 


Radius parameters (fm) 










Adopted value 


nucleus 


/ 2\l/2 

)ch '■0 




S'j2i%) 


Sh2i%) 


51/2 (%) 


5c2i%) 



T, = -1: 



10^ 


2.47(6) 


0.931(66) 


0.132(10) 


0.167(12) 


0.169(11) 


0.167(12) 


0.170(15) 


14q 


2.74(4) 


1.244(32) 


0.217(11) 


0.270(12) 


0.267(13) 


0.267(13) 


0.270(15) 


i«Ne 


3.00(3) 


1.361(20) 


0.251(6) 


0.386(9) 


0.387(8) 


0.381(10) 


0.390(10) 


^^Mg 


3.05(4) 


1.281(26) 


0.207(8) 


0.249(9) 


0.261(10) 


0.250 (8) 


0.255(10) 


28Si 


3.10(3) 


1.206(18) 


0.223(7) 


0.332(10) 


0.327(11) 


0.323(10) 


0.330(10) 


30 g 


3.24(2) 


1.223(13) 


0.812(15) 


0.728(15) 


0.730(17) 


0.750(16) 


0.740(20) 


34 Ar 


3.33(3) 


1.253(17) 


0.351(15) 


0.650(21) 


0.610(26) 


0.556(19) 


0.610(40) 


3«Ca 


3.48(2) 


1.269(10) 


0.402(11) 


0.727(17) 


0.674(18) 


0.596(12) 


0.710(50) 


«Ti 


3.60(5) 


1.316(22) 


0.359(14) 


0.563(26) 


0.572(29) 


0.578(33) 


0.555(40) 



T, = 0: 



26 m 


3.04(2) 


1.194(12) 


0.156(3) 


0.231(5) 


0.227(5) 


0.225(4) 


0.230(10) 


3*C1 


3.39(2) 


1.303(11) 


0.312(8) 


0.557(11) 


0.536(15) 


0.479(11) 


0.530(30) 


38m 


3.41(4) 


1.245(21) 


0.299(18) 


0.540(28) 


0.495(30) 


0.445(20) 


0.520(40) 


42Sc 


3.53(5) 


1.301(22) 


0.278(11) 


0.435(20) 


0.438(26) 


0.446(28) 


0.430(30) 


46y 


3.60(7) 


1.285(31) 


0.273(17) 


0.344(21) 


0.341(22) 


0.322(18) 


0.330(25) 


'^OMn 


3.68(7) 


1.260(30) 


0.315(20) 


0.439(27) 


0.455(33) 


0.438(28) 


0.450(30) 


54Co 


3.83(7) 


1.275(29) 


0.376(22) 


0.578(34) 


0.577(39) 


0.563(35) 


0.570(40) 


«2Ga 


3.94(10) 


1.271(42) 


1.31(11) 


1.10(11) 


1.07(11) 


1.01(8) 


1.05(15) 


«6As 


4.02(10) 


1.264(41) 


1.32(12) 


1.25(12) 


1.18(14) 


1.07(8) 


1.15(15) 


™Br 


4.10(10) 


1.264(39) 


1.43(13) 


1.11(13) 


1.03(14) 


0.85(6) 


1.00(20) 




4.18(10) 


1.276(37) 


0.68(9) 


1.51(14) 


1.38(18) 


1.20(12) 


1.30(40) 



of the wave function exactly matched that required for 
the known separation energy, Sp. With the well depth 
so fixed, we computed the radial wave functions for all 
proton states bound in that potential and constructed 
the charge density of the nucleus from the square of these 
functions: 



(23) 



nlj 



where 2j + 1 is the occupancy of protons in each orbital, 
nlj, and the sum is over the occupied orbitals. Here 



\R„ijir)\'r^dr/ / |i?„i,(r)|^r^dr, (24) 
Jo 



with Rnij (r) being the radial wave function of the proton 
with quantum numbers, nlj. We then determined the 
radius parameter of the Saxon- Woods potential, rg, by 
requiring the {r'^)lm computed from Eq. ( p3|) to match 
that determined from experimental electron scattering, 
Eq. (^. The value of tq is also given in Table and its 

error reflects the assigned error on {r'^)l.(j^ ■ 

In the shell model, the ^-particle wave functions, \JiTi) 
and \ JfTf), can be expanded into products of {A — 1)- 
particle wave functions |7r) and single-particle functions 



\j). In terms of this expansion, the Fermi matrix clement 



{T}MTfll\T,MT,\ 



Y^U{l\T,T^:\Tf)S''\^{\^■J)S''\f{\n■JW^ 



dr. 



(25) 



The expansion coefficients S^/^{i{\'K; j) and S^/^{f{\Tr;j) 
are generalised fractional parentage coefficients and rep- 
resent the spectroscopic overlap of the A- and {A — 1)- 
particle wave functions. The sum in Eq. (^5|) is over 
all parent states, |7r), and all single-particle orbitals ac- 
tive in the shell-model calculation. Note that the radial 
integrals, fij, are labelled with tt. These integrals are 
evaluated with eigenfunctions of the Saxon- Woods po- 
tential whose well depth is continually adjusted to match 
the separation energy to that particular parent state. If 
we do not allow the proton and neutron radial functions, 
i?P(r) and i?"(r), to vary with the parent states but fix 
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their asymptotic forms for all j to the separation energy 
of the ground state of the parent nucleus, then the sums 
over TT can be done analytically and the computed value 
of Sc'2 becomes independent of the shell-model effective 
interaction. Results of this calculation are given in Ta- 
ble and labelled ■ Results without this simplifying 
assumption are also given and labelled (5p^2- These lat- 
ter results depend on the effective interaction but not 
strongly. One reason for this is that in implementing 
Eq. (pa), we use experimental excitation energies in the 
(yl — 1) nucleus for the lowest-energy state of each spin 
and parity. The shell model is used to provide spectro- 
scopic amplitudes and the excitation energies of states in 
the (^—1) nucleus relative to the lowest state of that spin 
and parity. The difference between 6^2 indicates 

the role of the parentage expansions. 

So far, the two surface terms in Eq. ( pO| ) have not been 
included, Vg = 0, Vh = 0. It can be argued that the cen- 
tral part of the potential, which in principle should be de- 
termined from some Hartree-Fock procedure, should not 
be continually adjusted. Rather, any alteration should 
be to the surface part of the potential. Thus, in this 
method, we fix Vq separately for protons and neutrons 
to match the ground-state parent separation energies, Sp 
and Sn- For the excited parent states of excitation en- 
ergy. Ex, we adjust the strength of the surface term, Vg 
(keeping Vh = 0) so that the asymptotic forms match the 
separation energies Sp + and 5„ + Ex ■ These results 
are listed in Table ^ as 5^J2 ■ 

The second surface term, /i(r), is even more strongly 
peaked in the surface than g{r). Thus our fourth method 
is the same as the third, except that it is the second 
surface term, Vh, that was adjusted, keeping Vg = Q. 
These results are listed in Table ^ as (5^^. 

On average, the method III values of 5c2 are about 2% 
lower than the method II values; and method IV values 
about 7% lower than the method II values. These are 
not big differences. The errors on each individual en- 
try of 5c2 in the Table ^ reflects only the error in this 
quantity due to the uncertainty in the rms charge ra- 
dius (r^)^/^. Once again, as we have done in previous 
tables, the values tabulated for 5^2: ^C2^ ^C2 ^C2 
give the results from one particular calculation for each 
parent nucleus. Our adopted 6c2 values result from our 
assessment of all multiple-parentage calculations made 
for each decay - not just those shown in the preceeding 
three columns. The error on our adopted value reflects 
not only the uncertainty in the rms charge radius, but 
also the spread of results obtained with different shell- 
model effective interactions and the different procedures, 
II, III, and IV. 



D. Collected structure-dependent corrections: 
their reliability 

Our adopted values for the three nuclear-structure- 
dependent corrections, 5ns, 5ci and 5c2 are collected 



TABLE VI: Adopted values for the three nuclear-structure 
dependent corrections for superallowed Fermi /3 decay. 



Parent Sns{%) 5ci(%) 5c2(%) Sc ~ 5ms{%) 
T, = -1: 

10^ 
14q 

i«Ne 
22 Mg 
26Si 

30c 



*Ar 
^Ca 



= 0: 

26m 
34^1 

42Sc 

46 



50 



Mn 
'Co 
^Ga 
=As 
^Br 
'Rb 



-0.360(35) 
-0.250(50) 
-0.290(35) 
-0.240(20) 
-0.230(20) 
-0.190(15) 
-0.185(15) 
-0.180(15) 
-0.240(20) 



0.009(20) 
-0.085(15) 
-0.100(15) 

0.030(20) 
-0.040(7) 
-0.042(7) 
-0.029(7) 
-0.040(20) 
-0.050(20) 
-0.060(20) 
-0.065(20) 



0.010(10) 
0.050(20) 
0.230(30) 
0.010(10) 
0.040(10) 
0.195(30) 
0.030(10) 
0.020(10) 
0.220(100) 



0.040(10) 
0.105(20) 
0.100(20) 
0.060(30) 
0.095(20) 
0.055(20) 
0.040(15) 
0.330(40) 
0.250(40) 
0.350(40) 
0.130(60) 



0.170(15) 
0.270(15) 
0.390(10) 
0.255(10) 
0.330(10) 
0.740(20) 
0.610(40) 
0.710(50) 
0.555(40) 



0.230(10) 

0.530(30) 

0.520(40) 

0.430(30) 

0.330(25) 

0.450(30) 

0.570(40) 

1.05(15) 

1.15(15) 

1.00(20) 

1.30(40) 



0.540(39) 
0.570(56) 
0.910(47) 
0.505(24) 
0.600(24) 
1.125(39) 
0.825(44) 
0.910(53) 
1.015(110) 



0.261(24) 

0.720(39) 

0.720(47) 

0.460(47) 

0.465(33) 

0.547(37) 

0.639(43) 

1.42(16) 

1.45(16) 

1.41(21) 

1.50(41) 



in Table VI. Since their impact on the ft values is 
in the combination (Sc — S^s), (see Eq. (^)), where 
= Sci + 5c2, we list our results for this combination 
with the individual errors added in quadrature. Note 
that in the combination (Sc — Sns) all three corrections 
are in phase with the exception of the small Sns values in 
the cases of '^^'^Al and *^Sc. For the nine nuclei for which 
precision ft values have been measured, and ^''O of 
the -1 series, and ^e^Al to ^'^Co of the = se- 

ries, the nuclear-structure correction ranges from a low of 
0.26% for Al to a high of 0.72% for 38™K. Of partic- 
ular interest is that larger values are found at the upper 
end of the s,d-shell in the = —1 series and at the 
upper end of the p, /-shell in the — series. This is 
mainly due to the radial overlap correction, Sc2, which 
yields larger numerical values whenever a single-particle 
orbital with a radial node contributes importantly in the 
parentage expansions, such as the 2si/2 orbital in the up- 
per s,(i-shell and the 2^3/2, '^Pi/2 orbitals in the upper 
p, /-shell. 

There have been a number of previous calculations of 
Sc but only one of Sns- The latter was performed by 
one of the present authors using the same techniques 
described here but applied only to the nine well-known 
superallowed transitions and with similar - though dif- 
ferent in detail - shell-model calculations to ours; the re- 
sults for those transitions are very similar to the present 
results, well within the error bars in all cases. 
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The more numerou s re sults from previous 5c calcula- 
tions appear in Table VII, where they are compared with 
our present results. Four groups of authors have pub- 
lished values for 5c, the first in 1973. In the table, we 
present the most recent results from each group for each 
transition. The values in the first column are those cal- 
culated previously by us, reported first in refs. |2q] 
and then refined in more recent publications \^ p6[ . 
These were based on the same methods as those used 
here: shell-model calculations to determine 5ci, and full- 
parentage expansions in terms of Woods-Saxon radial 
wave functions to obtain 5c2- Ormand and Brown, whose 
values for 5c appear in column 2, also employed the 
shell model for calculating Jci, but they derived 5c2 from 
a self-consistent Hartree-Fock calculation. Both of these 
independent calculations for 5c ~ those in columns one 
and two - reproduce the measured coefficients of the rele- 
vant isobaric multiplet mass equation, the known proton 
and neutron separation energies, and the measured ft- 
values of the weak non-analogue 0+ 0+ transitions ||] 
where they are known. The agreement of these calcula- 
tions with our new results is rather good, especially for 
the well known nine. In the cases of the less well known 
Tz = —1 nuclei between ^*Ne and ^^Ti, the differences 
are in general larger, but this reflects improvements to 
sc?-shell calculations realized since 1973, when the only 
previous calculations j25| were published. 

The other two previous calculations shown in the table 
provide a valuable check that these 5c values do not suf- 
fer from severe systematic effects. Sagawa, van Giai and 
Suzuki 1^ have added RPA correlations to a Hartree- 
Fock calculation that incorporates charge-symmetry and 
charge-independence breaking forces in the mean-field 
potential to take account of isospin impurity in the core; 
the correlations, in essence, introduce a coupling to the 
isovector monopole giant resonance. The calculation is 
not constrained, however, to reproduce known separation 
energies. In addition, the authors themselves|Q admit 
that their HF+RPA calculations cannot properly take 
account of pairing in open-shell nuclei; as a consequence, 
the discrepancies between their values and the others for 
■^'^Cl and ^Smj^ jg j^q^ considered significant. Clearly the 
overall trend of the shcU-model-based calculations is well 
reproduced by these very-different calculations, thus rul- 
ing out the possibility that the former had missed sig- 
nificant core contributions. Finally, a large shell-model 
calculation has been mounted for the A = IQ case by 
Navratil, Barrett and Ormand |2^. This "microscopic" 
calculation of 5c also supports the results of the more 
macroscopic calculations reported here and in columns 1 
and 2. 

We can now address the question of whether the CKM 
unitarity problem might be removed by plausible changes 
in the calculated structure-dependent corrections embod- 
ied in 5c — 5ms- As can be seen from Table VII, the 
typical value of 5c —5ns is of order 0.5% for the nine 
well-known cases currently used in the unitarity test. To 
remove the unitarity problem, the nuclear-structure de- 



TABLE VII: Calculated values for the isospin symmetry- 
breaking correction, 5c in percent units. Previous calcula- 
tions are compared with the present results. 



Parent 


Towner 


Ormand 


Sagawa 


Navratil 


Present 




o6 LLcLL try 


06 JDiOWli 


CL U/i 


Ci (It 


work 


= -1: 














0.18(2) 


0.15(9) 


0.00 


0.12 


0.18(2) 


14q 


0.28(3)" 


0.15(9) 


0.29 




0.32(3) 


i«Ne 


0.45(3) 








0.62(3) 




0.35(3) 








0.27(2) 


28 Si 


0.42(4) 








0.37(2) 


30g 


1.21(10) 








0.94(4) 


34 Ar 


1.04(9) 








0.64(4) 


3^Ca 


0.89(9) 








0.73(5) 




0.62(6) 








0.78(11) 


= 0: 














0.33(5)" 


0.30(9) 


0.27 




0.27(2) 


3*C1 


0.64(7)" 


0.57(9) 


0.33 




0.64(4) 


38jTi 


0.64(7^' 


0.59(9) 


0.33 




0.62(5) 




0.40(6)^ 


0.42(9) 


0.44 




0.49(4) 


46y 


0.45(6/ 


0.38(9) 






0.43(3) 


50Mn 


0.47(9/ 


0.35(9) 






0.51(4) 


s^Co 


0.61(6/ 


0.44(9) 


0.49 




0.61(4) 


^^Ga 




1.26-1.32S 


1.42 




1.38(16) 


""As 




1.41-1.63» 


0.78 




1.40(16) 


™Br 




1.11-1.41" 






1.35(21) 


^4Rb 




0.91-1.05« 


0.74 




1.43(40) 



"Both Sci and 5c2 are taken from Towner, Hardy and Harvey||9|], 
except as noted. 
''Both Sci and Sc2 are taken from Ormand and Brown 
'^SGII results from Sagawa, van Giai and Suzu ki [pq ] 
''Value of Sq from Navratil, Barrett, Ormand pol 
■^The values of 5ci are taken from ref. [26|| 
•^The iaj ues of Sci are taken from ref. 

*Ref. pvj uses two methods to calculate Sc2 for these cases; to be 
consistent with other numbers in this column, we quote the results 
for Hartree-Fock wave functions. 



pendent corrections, {5c — 5ns): would all have to be 
raised to around 0.8%. Neither the present work nor 
any previous calculation gives any indication that such 
a systematic shift is plausible under any reasonable cir- 
cumstances. 

The structure-dependent corrections have another 
more impressive credential, one that is not often appreci- 
ated: they are demonstrably effective in bringing the dis- 
parate experimental /t- values into agreement with CVC. 
If the experimental /t-values were left uncorrected, their 
scatter would be quite inconsistent with a single value 
for the vector coupfing constant, Gv. Once corrected, 
the resulting JFt-values are in excellent agreement with 
this expectation. In a very real sense, it can be said that 
CVC supports the structure-dependent corrections. This 
point will be amplified in the next section. 
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TABLE VIII: Calculated values for the corrected Tt-vahies 
based on the adopted (average) Sc values and world-average 
experimental ft values. 



Parent 


Adopted 


Experimental 


Corrected 


nucleus 


oc(/o) 






T, = -1: 










0.17(3) 


3038.7(45) 


3072.7(48) 


14 Q 


0.24(3) 


3038.1(18) 


3069.4(26) 


Tz = 0: 










0.29(3) 


3035.8(17) 


3071.4(22) 


3*01 


0.61(3) 


3048.4(19) 


3070.6(25) 


38m 


0.61(3) 


3049.5(21) 


3070.9(27) 


^^Sc 


0.46(3) 


3045.1(14) 


3075.7(24) 


46y 


0.41(3) 


3044.6(18) 


3074.4(27) 




0.43(3) 


3043.7(16) 


3072.9(28) 


54 Co 


0.53(3) 


3045.8(11) 


3072.1(27) 






Average Tt 


3072.2(8) 








0.6 



"Average of present results with those_aLOrmand and Brown |p^; 
both are listed individually in Table VII. The uncertainties are 
explained in the text. 

''Data taken from ref. Illl. 



graphically in Fig. |^. The upper panel shows the uncor- 
rected experimental ft values and the lower panel the 
corrected !Ft values with the average indicated by a hori- 
zontal line. Evidently, in these cases, at the current level 
of precision the nucleus-dependent corrections act very 
well to remove the considerable "scatter" that is appar- 
ent in the experimental /i- values and is effectively absent 
from the corrected ^i-values. As mentioned already, the 
consistency of the corrected Tt values {x^/i^ — 0.6) is a 
powerful validation of the calculated corrections used in 
their derivation. 

Of course it is only the relative values of {6c—Sns) that 
are confirmed by the absence of transition-to-transition 
variations in the corrected ^t-values. However, 6c itself 
represents a difference - the difference between the par- 
ent and daughter-state wave functions caused by charge- 
dependent mixing. Thus, the experimentally determined 
variations in Sc are actually second differences. It would 
be a pathological fault indeed that could calculate in de- 
tail these variations {i.e. second differences) in Sc while 
failing to obtain their absolute values (i.e. first differ- 
ences) to comparable precision. 

We have argued that decreasing the radiative correc- 
tion S'j^ from 1.4% to 1.1%, or Al from 2.4% to 2.1% is 



III. THE J'^t- VALUES: PRESENT STATUS AND 
FUTURE PROSPECTS 

With improved calculations for Sc and Sns, we are 
now in a position to extract corrected J^t- values from the 
current world data for the nine well known experimental 
/t-values. To do so, we follow the same procedure we 
have used in the past[0, |2) to arrive at values for Sc that 
best represent the results from the two groups that have 
made complete calculations: in the present situation that 



means we use Table VII and take an unweighted averaa 



of the results in column three (Ormand and Brown p7|) 
with those in column six (present work). Noting that 
there is a small systematic difference of 0.08% between 
the two sets of calculations, we remove that difference 
and then analyze the scatter of all nine pairs of Sc re- 
sults about their respective averages to obtain a standard 
deviation of 0.034%. Our adopt ed Sc values appear in 
the second column of Table VIIl where they also include 
the adopted "statistical" uncertainty of 0.034%. (The 
"systematic" uncertainty of ±0.04%, obtained from the 
average difference between the two calculations of Sc, 
need not be applied until Gv is extracted from the J-t- 
values.) 



The next columns in Table VIII contain the experimen- 
tal /i- values, which we have simply taken from ref. [Q, 
and the corrected J-"t-values, which we have calculated 
from Eq.(^) using Sc from the first column of this table, 
Sns from column two of Table VI and S'j^ from the last 
column of Table |. The average JFi- value and the corre- 
sponding P6r degree of freedom also appears at the 
bottom of the table. The same information is presented 



ft 



3052 
3050 
3048 
3046 
3044 
3042 



'"C ^"cr "Sc ^"Mn 



il 



'If' 



i 



5 10 15 20 25 30 



7t 



3080 
3078 
3076 
3074 
3072 



C "CI Sc "Mn 

"O ""Al '""K "V "Co 



5 10 15 20 25 30 

Z of daughter 

FIG. 1: Comparison of experimental ft-vahies and the cor- 
rected .T^t-values for the nine well-known superallowed tran- 
sitions. This illustrates the effect of the calculated nucleus- 
dependent corrections, which change from transition to tran- 
sition. (The effect of S'^ is virtually the same for all cases.) 
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unlikely to be the solution to the CKM unitarity prob- 
lem; and that there is no support from calculations for an 
average increase in the nuclear-structure dependent cor- 
rection, {Sc — Sns), from 0.5% to 0.8%. We are therefore 
confident that the unitarity result in Eq. , which is un- 
changed by our new calculations, incorporates structure- 
dependent corrections that are correct within their stated 
uncertainties. Nevertheless, these uncertainties are con- 
servatively assigned and, as we remarked in the intro- 
duction, they contribute significantly to the overall un- 
certainty of the unitarity test. There is every reason to 
continue to focus on these corrections, both experimen- 
tally and theoretically, with a view to reducing their un- 
certainties still farther. 

One way to do so, of course, would be to increase the 
precision of the /t-values for the nine cases tabulated in 
Table VII] and thus improve the comparison with CVC 
that is illustrated in Fig. ^ However, given the large 
amount of high-quality data that is already incorporated 
in these nine /i-values, significant improvements are un- 
likely in the near term. A more promising experimen- 
tal approach to testing Sc is offered by the possibility 
of increasing the number of superallowed emitters ac- 
cessible to precision studies. Two series of 0+ nuclei 
present themselves: the even-Z, = —1 nuclei with 
18 < A < 42, and the odd-Z, = nuclei with A > 62. 
The main attraction of these new regions is that the cal- 
culated values of Sc — Sns for the superallowed transi- 
tions are larger, or show larger variations from nuclide 
to nuclide, than the values applied to the nine currently 
well-known transitions (see Table VI). In principle, then, 
they afford a valuable test of the accuracy of the Sc cal- 
culations. It is argued that if the calculations reproduce 
the experimentally observed variations where they are 
large, then that must surely verify their reliability for 
the original nine transitions whose Sc values are consid- 
erably smaller. The calculations reported here, the only 
complete set available for all these new cases, should pro- 
vide a sound basis to which new experimental data can 
be compared. 

Currently, the greatest attention is being paid to the 
Tz = emitters with A > 62, since these nuclei are being 
produced at new radioactive-beam facilities, and their 
calculated Sc corrections had previously been predicted 
to be large 



7, 28 



It is likely, though, that the required 
experimental precision will take some time to achieve. 
The decays of these nuclei are of higher energy and each 
therefore involves numerous weak Gamow- Teller tran- 
sitions in addition to the superallowed transition [p9[ . 
Branching-ratio measurements will thus be very demand- 
ing, particularly with the limited intensities likely to be 
available initially for most of these rather exotic nuclei. 
In addition, their half-lives are considerably shorter than 
those of the lighter superallowed emitters; high- precision 
mass measurements (±2 keV) for such short-lived activ- 
ities will also be very challenging. 

More accessible in the short term will be the = — 1 
superallowed emitters with 18 < A < 42. There is 



good reason to explore them. For example, the calcu- 
lated value of {Sc — Sns) for '^"S decay, though smaller 
than those expected for the heavier nuclei, is actually 
1.13% - larger than for any other case currently known - 
while ^^Mg has a low value of 0.51%. If such large differ- 
ences are confirmed by the measured /t-values, then it 
will do much to increase our confidence in the calculated 
Coulomb corrections. To be sure, these decays will also 
provide a challenge, particularly in the measurement of 
their branching ratios, but the required precision should 
be achievable with isotope-separated beams that are cur- 
rently available. 



IV. CONCLUSIONS 

We have presented a new and consistent set of calcula- 
tions for the nuclear-structure dependent corrections, {Sc 
- ^Ns)i required in the analysis of superallowed 0"*" — > 0^ 
beta decay. Twenty transitions have been included in our 
calculations, the nine well known ones already used in the 
CKM unitarity test, and eleven more that are likely to 
be accessible to precise measurements in the future. The 
unitarity test itself is unchanged by our calculations, one 
of several indications we offer that these corrections are 
under control within their stated uncertainties. We have 
also argued that the structure-independent radiative cor- 
rections are similarly sound. If the apparent deviation 
from unitarity is to be resolved without demanding some 
extension to the Standard Model, the only remaining pos- 
sibility is through undiscovered errors in Vus , whose value 
is currently derived from Kg^ decay [Q, ^ and has not 
been revisited in nearly 20 years. 

We have also shown that the uncertainty quoted for the 
unitarity test can most effectively be improved by reduc- 
tions in the uncertainties of and {Sc — Sns)- We have 
outlined an experimental method by which the latter can 
be improved, and have provided the full set of calculated 
corrections that can be tested against experiment. The 
stage is now set for a new infiux of experimental results 
on previously unexplored superallowed transitions, from 
which the calculated structure-dependent corrections can 
be tested and confirmed or refined. In cither case, the 
uncertainties should be reduced and the unitarity test 
sharpened. 
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APPENDIX A: EFFECTIVE INTERACTIONS 

In the tables of results presented in the main text, 
we have only provided one set of values for each decay 
studied. However, for each nucleus, many calculations 
were performed with varying choices of effective interac- 
tions and shell-model spaces. The error assigned to the 
adopted values reflects the spread in the results and our 
estimate of the uncertainty in the calculated value based 
on the quality of the shell-model calculation. 

The choice of an effective interaction is easily made 
for shell-model calculations in light nuclei whose princi- 
pal configurations involve several valence nucleons away 
from major shell closures. There are well established in- 
teractions that give excellent fits to spectra. For A = 10, 
we use the Cohen-Kurath||3l[ interaction, (8-16)POT, 
and for A = 22, 26, 30 and 34, we use the univer- 
sal s, c?-interaction, USD, of Wildenthal|3^]. For nuclei 
with A = 46, 50 and 54, we considered two interactions: 
the Kuo-Brown G-matrix|Q as modified by Poves and 
Zuker[Q and denoted KB3, and the /p- model indepen- 
dent interaction of Richter et aZ.js^ and denoted FPMI3. 
For nuclei with A — 50 and 54 it was not possible to per- 
form untruncated calculations in the full fp space; our 
calculations only contain if7/2T~''iP3/2, f5/2,Pi/2T con- 
figurations with r < 2. In this truncated calculation, the 
spectrum obtained for 0+ states in A = 50 and 54 is in 
very poor agreement with experiment, a much larger en- 
ergy gap between the ground state and first excited 0^ 
being obtained. Thus, we have made further adjustments 
to the interaction centroids to obtain a much improved 
spectrum in the truncated space. 

For nuclei with A = 62, 66 and 70 we considered the 
model space (p3/2j /5/2)?'i/2)"i with n = ^ — 56, which is 
based on a closed /7/2 shell at the ^^Ni core. This model 
space is the one used by Koops and Glaudemans [|6) in 
their study of nickel and copper isotopes. We found this 
model space, with a modified surface delta interaction 
(MSDI) as used in j36|, gave acceptable spectra for the 
beta-decaying nuclei, with excited 0^ states at about the 
right excitation energy. 

The problem cases were A = 14, 18, 38, 42 and 74. In 
each of these cases, the experimental excited O'^ states 
are at a much lower energy than can be obtained in shell- 
model calculations. This is symptomatic of the presence 
of deformed configurations intruding among the spherical 
shell-model configurations. For example, in the A = A2 
spectrum the lowest-energy states are predominantly two 
particles outside a closed ^°Ca core, \2p), but lying low 
in the spectrum are 'intruder' states with a configuration 
of four particles and two holes, |4p-2/i). Mixing between 
these configurations must occur, and it is difficult to ob- 
tain the correct degree of mixing with the shell model. 
Shell- model calculations that attempt to mix \2p) and 
4p-2/i) configurations encounter what has been called 



p7[ the "nhuj catastrophe". The presence of |4p-2/i) con- 
figurations depresses the \2p) states, opening up a large 
energy gap between the \2p) and |4p-2/i) states. This 
would be corrected somewhat if the model calculation in- 
cluded |6p-4/i) states as well, since the role of the \6p-4:h) 
states is to depress the |4p-2/i) states. Thus if the model 
space is truncated to include only \2p) and |4p-2ft,) states, 
the depression driven by the \6p-4h) states on the |4p-2/i} 
states is absent. In an attempt to circumvent this catas- 
trophe we weakened the cross-shell interactions. Specif- 
ically, at mass 14, 18, 38 and 42 we used the Millener- 
Kurath |3^ interaction to evaluate the {2p\V\Ap-2h) ma- 
trix elements. We multiplied these matrix elements by 
a factor, /, that ranges from 0.0 to 1.0. When / = 0.0, 
there is no mixing between \2p) and \Ap-2h) configura- 
tions, and when / = 0.6 the ground-state wave function 
is approximately 80%|2p} and 20%\Ap-2h). Our strategy 
was to adjust / so that the excited 0+ energy is approx- 
imately equal to the experimental excitation energy. We 
have examined the sensitivity of our results to variations 
in / and ensured that the spread of values obtained were 
within the assigned errors attributed. 

There are some older interactions that operate in very 
restrictive model spaces, but remove the nhto catastro- 
phe by allowing mixing between \2p), |4p-2/i) and \6p-Ah) 
configurations. These are the Zuker-Buck-McGrory|p9| 
interaction, ZBM, as modified by Zuker [^, which uses 
the Pi/2, Si/2 and d^/2 orbitals for the A = 14 and 18 nu- 
clei; and the Federman-Pittel[^ interaction, FP, which 
uses the ^3/2 and fj/2 orbitals for the A = 38 and 42 
nuclei. 

Finally, at mass 74 there is a related but slightly dif- 
ferent problem. The spectrum in a (^3/25/5/2:^1/2)^^ 
model space gives about the right density of natural- 
parity states. The difficulty is the presence of unnatural- 
parity states lying low in the spectrum (for example, ^^Br 
has a 5/2+ at only 280 keV excitation, while ^^Rb has a 
probable 3/2+ at 40 keV). Further, the excited 0+ state 
in '^'^Ki is at only 508 keV, whereas the (P3/2, /5/2,Pi/2)^^ 
model calculation puts the state at 2550 keV. The infiu- 
ence of the I59/2, 2^5/2 and possibly 137/2 orbitals is 
evidently quite strong at the end of the p, f shell. Thus, 
we have used the following model space 



(P3/2, /5/2,Pl/2 



,18 



- (P3/2)*(/5/2,Pl/2)^(59/2, 0^5/2 , 57/2)^ ■ 

(Al) 



Let us call the first term in Eq. (Al) the OHlo term, and 
the second term with two nucleons promoted to the d, g 
shell the 2huj term. Because of the "n?ia; catastrophe", 
we again multiply all {0hio\V\2hLu) matrix elements by 
a factor, /, and adjust / so that the excited 0+ state 
in ^''Kr is reproduced at its experimental location. All 
matrix elements were then calculated with the MSDI 
interaction[06|. 
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